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LOOP STATISTICS IN THE TOROIDAL HONEYCOMB DIMER 

MODEL 

By Cedric Boutillier and Beatrice de Tiliere 1 
UPMC University Paris 06 and CNRS and Universite de Neuchdtel 

The dimer model on a graph embedded in the torus can be inter- 
preted as a collection of random self-avoiding loops. In this paper, 
we consider the uniform toroidal honeycomb dimer model. We prove 
that when the mesh of the graph tends to zero and the aspect of 
the torus is fixed, the winding number of the collection of loops con- 
verges in law to a two-dimensional discrete Gaussian distribution. 
This is known to physicists in more generality from their analysis 
of toroidal two-dimensional critical loop models and their mapping 
to the massless free field on the torus. This paper contains the first 
mathematical proof of this more general physics result in the specific 
case of the loop model induced by a toroidal dimer model. 

1. Introduction. Two-dimensional critical loop models are believed by 
physicists to renormalize at criticality onto a Gaussian free field theory 
(Coulomb gas). On this basis, they are able to derive explicit formulae for 
partition functions on the torus [1, 3, 6, 7, 15], from which they obtain infor- 
mation about the asymptotic distribution of the sum of the winding numbers 
of the loops. In this paper, we give the first mathematical proof of this result, 
in the case of the toroidal uniform dimer model on the honeycomb lattice. 

Let G = (V(G), E(G)) be a graph embedded in the torus. A dimer con- 
figuration or perfect matching M of G is a subset of edges of G such that 
every vertex of G is incident to exactly one edge of M. Superimposing M 
onto a reference dimer configuration Mq yields a collection of self-avoiding 
loops together with doubled edges; see Figure 2. Loops in this collection can 
wind around the torus horizontally and vertically. 

We consider the case where G is a quotient of the hexagonal lattice and 
where dimer configurations are chosen uniformly among all dimer config- 
urations of G. The algebraic sum of the winding numbers of the loops in 
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the collection is then a random variable. The main result of this paper can 
loosely be stated as follows (see Theorem 2 of Section 1.3 for a precise state- 
ment): 

Theorem 1. When the mesh of the graph G tends to and the aspect 
of the torus is fixed, the probability that the algebraic sum of the winding 
numbers of the loops has horizontal component k and vertical component I 
converges to 

_}_ ^n(k 2 /p+ P e 2 )/2 

where Z p = J2(k,e)ez 2 e _7r ( fc2//p+p ^ 2 )/ 2 is the normalizing factor and p is the 
ratio of the two side lengths of the torus. 

Let us now describe the setting in more detail. 

1.1. Toroidal honeycomb dimer model. The dimer model is a statisti- 
cal mechanics model introduced to represent the adsorption of diatomic 
molecules on the surface of a crystal. When, in addition, the underlying 
graph is bipartite (as is the case for the square lattice or the honeycomb lat- 
tice), this model can be interpreted as a random interface model in dimension 
2 + 1, via the height function [17]. The dimer model has the attractive fea- 
ture of being exactly solvable [8, 9, 16] and is believed to be conformally 
invariant in the scaling limit; for rigorous results on this, see [5, 10, 11, 12]. 

In the case where the surface of the crystal is modeled by the regular 
hexagonal lattice H, we speak of the honeycomb dimer model. The honey- 
comb lattice H has a natural embedding in the plane, in which all faces 
are regular hexagons of side length 1 (or, equivalently, the dual faces are 
equilateral triangles of side length v3). Consider the two vectors x and y 
represented in Figure 1. The lattice H and its bipartite coloring are invari- 
ant under the action of x and y by translation. If, for every m,n € N*, we 
define L mjn to be the group of translations spanned by mx and ny, then 
l^ m ,n is a subgroup of the symmetry group of H . The toroidal graph H m ^ n 
is defined to be the quotient H m>n = HfL m>n . 

The graph H m ^ n can be obtained by cutting out an mx x my rectangle in 
H and then gluing opposite sides together. A horizontal side of the rectangle, 
oriented from left to right, is mapped to an oriented closed loop on the torus, 
as is a vertical side of the rectangle, oriented from bottom to top. Let us 
denote by j h and j v these two oriented closed loops. The first homology 
group H\(T) is spanned by the homology classes of j h and j v . 

The modulus of the torus is the ratio of the two complex numbers rep- 
resenting the vertical and the horizontal sides of the rectangle. Since the 
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length of the horizontal side is 3m and the vertical side has length nv3 by 



construction, the modulus of H mn is i-j=— . 

A dimer configuration of H m>n is a perfect matching of H mn , that is, a 
subset of edges M of H m ^ n such that every vertex of H m>n is incident to 
exactly one edge of M. Let us denote by J\4(H mn ) the set of dimer con- 
figurations of H m>n . Suppose that a positive weight function v is assigned 
to edges of H mjn , that is, every edge e has weight v{e). Then, every dimer 
configuration M of iT m , n has an energy S(M) = — 2~Z e eMl°S^( e )- The prob- 
ability of occurrence of the dimer configuration M of H m n is given by the 
Boltzmann measure [x m ,n- 

a m-— - n ^ (e) 

where Z m ,„(i/) = EMeAl(iJ ro ,„) U e eM u ( e ) is tne normalizing constant, known 
as the partition function. 

When the weight function v is periodic, it is known, using classical subad- 
ditivity arguments (as in [4]) that the quantity —^logZ niin (u) converges 
when m and n tend to oo. The limit is denoted by f(z^) and is called the free 
energy per fundamental domain. 

In this paper, we consider dimer configurations of H m ^ n chosen with re- 
spect to the uniform measure (i.e., the Boltzmann measure corresponding 
to weights 1 on all edges). The corresponding partition function is called 
the uniform partition function. In this case, the free energy per fundamental 
domain, simply denoted by f, is given by the following formula [8], whose 
derivation is recalled in Section 3.2: 

(1) f = "A / 2 "log(2(cos^ + l)-e^)#dV. 

4tt z Jo Jo 

1.2. Toroidal dimer model and self-avoiding loops. The dimer model on 
the toroidal graph H m n can be interpreted as a collection of self-avoiding 




Fig. 1. The lattice H is invariant under the action of x and y by translation (left). 
Fundamental domain H\,\ (right). 



4 



C. BOUTILLIER AND B. DE TILIERE 




Fig. 2. The superimposition of M and Mo consists of doubled edges and alternating 
loops. 

loops as follows. Let Mq be a fixed dimer configuration of H mn and let 
M be any other dimer configuration of H m ^ n . The superimposition of Mq 
and M then consists of self-avoiding doubled edges and alternating loops, 
where doubled edges are edges covered by a dimer in both Mq and M, and 
alternating loops are cycles whose edges are, in alternation, dimers of Mq 
and M; see Figure 2. This feature is due to the fact that, by definition of 
perfect matchings, every vertex of H m ,n is incident to exactly one edge of 
Mq and one edge of M. 

Orienting the dimers of M from their white end to their black end and 
the dimers of Mq from their black end to their white end gives rise to an 
orientation of the loops. Let us denote by M Q Mq the set of oriented loops 
obtained from this superimposition; see Figure 3 for an example. 

A loop C of M Q Mq can then be seen as a closed path on the torus T. 
The equivalence class [C] of C in the first homology group H\ (T) ~ Z 2 can 




Fig. 3. The superimposition M Q Mo consists of oriented loops. In this example, the 
winding number is windjif (M) = (1,0). 
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Fig. 4. Dimer configuration of H\,z which generates the reference matching Mq. 



be decomposed in the basis ([y h ], [j v ])- Its coordinates (C h , C v ) in this basis 
are called the winding number of C: C h (resp., C v ) is the algebraic number 
of times the loop C winds horizontally (resp., vertically) around the torus. 
The winding number of the dimer configuration M, denoted by windM () (M), 
is the sum of the winding numbers of all of the loops contained in M Q Mq: 

wind Mo (M)= 

C loop 
in MQM 

An example of a computation of windjif (M) is given in Figure 3. Note that 
the dependence on Mq of the winding number of M is quite simple: if M\ 
is another dimer configuration, then 

wind M i {M) = wind Mo (M) - wind A / (Mi). 

From now on, we assume that n is a multiple of 3 and fix the reference 
dimer configuration Mq of H m ^ n to be the one generated by translations of 
the dimer configuration of H\ t z °f Figure 4. We also drop the subscript Mq 
in windMo(M). 

Since dimer configurations of H m ,n are chosen according to the uniform 
measure, for every m,n, wind(-) is a random variable. Let us call it the 
winding number and, in order to stress the dependence on m and n, denote 
it by wind mi „(-). 

The main result of this paper is an explicit expression for the asymptotic 
distribution of the random variables (wind m) n) when m,n tend to infinity 
and when the modulus i—^— of H m n converges to in for some p > 0. 

Vim ).■_.! 

1.3. Statement of result. Let us assume that the modulus ^"^r — of ^m^n 
converges to ip for some p > 0. Recall that when n is a multiple of 3, wind mjn 
is the winding number of the uniformly distributed dimer configurations 
of H mt1l , computed with respect to the reference dimer configuration Mq 
defined above. Then, the main result of this paper is the following theorem: 
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Theorem 2. In the joint limit m,n — > oo, — ► i/te sequence of 
random variables (wind m , jn ) converges in distribution to the two-dimensional 
discrete Gaussian random variable windp whose law is given by 

(2) V{k,£)£Z 2 F[wmd p = (k,£)} = ^e~ 7T{k2/p+p£2)/2 , 

where Z p = Z ik ,i)^e~ n{k2/p+pe2)/2 - 

• A similar result was obtained by Kenyon and Wilson [14] in the case of 
the square lattice embedded in the cylinder. Working on the torus makes 
computations much more difficult since it means dealing with the toroidal 
partition function in the proof, which is a combination of four terms [9] 
(instead of one, as in the cylinder case). Moreover, we have to extract 
information about the two components of the winding number (instead of 
one, as in the cylinder case). Note, also, that in proving Theorem 2, we 
give a full asymptotic expansion of a perturbation of the uniform partition 
function; see Theorem 4, Section 2 below. 

• It was brought to our attention, after the acceptance of this paper, that 
an asymptotic expansion of the uniform partition function was obtained 
in the physics literature by Ferdinand, in the case of the square lattice [7]. 
Nevertheless, let us stress the following facts: the expansion of [7] is not 
perturbative, it is not done to the same level of mathematical rigor and 
no information about the distribution of the winding number is inferred. 

• Theorem 2 holds if the honeycomb lattice is replaced by the square lat- 
tice. We believe that the techniques and ideas applied here could be used 
to extend the result to quotients of the honeycomb and square lattice, 
embedded on tori with modulus in the upper half complex plane (and not 
just ip, p > 0). More generally, we conjecture the result to be true when 
considering the dimer model on any periodic bipartite graph within the 
liquid phase [13], with an appropriate embedding. 

1.4. Outline of the paper. 

• In Section 2, we prove that the moment generating function of wind mjn 
can be expressed in terms of a perturbed uniform partition function 
Z m ,n{oi,0)\ see Section 2 for a definition. Theorem 4 then gives a full 
asymptotic expansion of Z m>n (a, /?), from which we deduce pointwise con- 
vergence of the moment generating function of wind. m)Tl to the moment 
generating function of windp and, hence, Theorem 2. 

• The remainder of the paper consists of the proof of Theorem 4, giving 
the full asymptotic expansion of the perturbed uniform partition function 
Z mn (a, /3). For the reader's convenience, the proof is split into two parts: 
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- By [8], the partition function Z m ^ n {a,(3) can be expressed as a linear 

combination of four terms, Zm$ (a, j3), a,n G {0,1}. Proposition 8 of 
Section 3 gives a full asymptotic expansion for each of the four terms 
as a function of Jacobi's four elliptic theta functions. The 
proof of Proposition 8 is postponed until Section 4. Proposition 9 then 
gives an explicit expression for the combination of Jacobi theta functions 
involved in the expression of Z mtn (a, (5). The proof of Theorem 4 is thus 
completed, apart from the proof of Proposition 8. 

- Section 4 consists of the proof of Proposition 8, giving the full asymp- 
totic expansions of the four terms Zm V n{a,j3) as functions of Jacobi's 
four elliptic theta functions. 

2. Winding number and partition function. Define F m ^ n to be the mo- 
ment generating function of the random variable wind mjn : 

V(q, 13) G R 2 F m>n (a, p) = E[e~™ indm <" < a ' /3) ] 

= ]T P[wind m , n = (M)]e- w(afc+/ ^. 

(M)6Z 2 

Convergence in distribution of the sequence (wind mj „) to the two-dimensional 
Gaussian random variable wind p given in (2) is equivalent to pointwise con- 
vergence of the sequence (F mtn ) to the corresponding moment generating 
function. 

Lemma 3 below gives an expression of F m ^ n in terms of a perturbed uni- 
form partition function, defined as follows. Let us introduce the appropriate 
choice of perturbed edge- weights. Define edges to be of type I (resp., II, III), 
as in Figure 5. Then, for a.,/3 G M, let us assign weights a = e - an /( 2m ) to 
edges of type I, 6 _1 = e~^^ 2n ^ to edges of type II and b = e^ 71 "^ 2 ™) to edges 
of type III. In the sequel, we shall also use the notation A = a~ 2m = e a7T , 
B = b 2n = . Observe that the weights a, b tend to 1 when m, n tend 
to infinity and thus yield a perturbation of the uniform partition function. 
They are used to collect information on the uniform measure. Let us denote 
by Z mjn (a,P) the partition function of the graph H miTl corresponding to 
these weights and let us call it the perturbed uniform partition function or, 
in short, perturbed partition function. 

Lemma 3. The moment generating function F m!n (a, (3) o/wind min and 
the perturbed uniform partition function Z m , n (a,f3) are related in the fol- 
lowing way: 

K mA a,/J)-e ZmAm - 
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Proof. For i =1, II, III, let Ni(M) be the number of edges of type i 
in the dimer configuration M . Then, by definition, the partition function 

Z m!n (a,f3) is 

Z mn (a,/3) = V ^ e -a7r/(2m)^A r i(M)^ e -/3 7r /(2n)- ) A r ii(M) ( - e /3 7 r/(2n)^ni(Af) 
MeM{H mtn ) 

= V A -N I (M)/(2m) B (N lu (M)-N u (M))/(2n) 
MeM(H mtn ) 

where, we recall, A = e 71 "", B = e 71 "^. 

Let us compute the two components windj^ n and windj^ n of wind m)n as 
functions of Ni,Nu,Nm. Consider the 2n left-to-right horizontal paths of 
the dual graph -ff^ru as m Figure 6 (left). Then, for any such path 7, 
windJ^ n (M) is equal to the number of positive (right-to-left) crossings of 
M Q Mq along 7 minus the number of negative (left-to-right) crossings of 
M Q Mq along 7. Summing over all 2n horizontal paths and observing that 
jV m (Mo) = (M ), we obtain 

2n • wind^ n (M) = N U (M) + N m (M ) - N m (M) - N U (M ) 

= N u (M)-N m (M). 

In a similar way, considering the 2m top-to-bottom vertical paths of Figure 
6 (right) and observing that Ni(Mq) = 2mn/3, we obtain 

2m • wind^ n (M) = Ni(M) - Ni(M ) = Nj{M) - ——. 




Fig. 5. The three types of edges around a white vertex (left) and around a black vertex 
(right). 




LOOP STATISTICS IN THE TOROIDAL HONEYCOMB DIMER MODEL 9 

Plugging this into the partition function Z mtn (a,/3), we obtain 

Z m , n (a,p) = A~ n / 3 ^ wind --( M ) J B- wind ™."( M ) 
MeM(H m , n ) 

= e -nan/3 C k ^^ ak+ ^\ 

where Ck,e is the number of dimer configurations whose winding number is 
(k,£). The proof is completed by recalling that dimer configurations of H m n 
are chosen with respect to the uniform measure, which implies that 

P[wind m , n = (M)] 



Z m , n (0,0) 

[Z m ^ n (0,0) is the uniform partition function]. 

□ 

Theorem 4 below gives a precise asymptotic expansion of Z m>n (a,P). 
Combined with Lemma 3, this yields pointwise convergence of the sequence 
{Fm,n) to the moment generating function F p of the discrete Gaussian ran- 
dom variable windp and, hence, Theorem 2. 

Theorem 4. In the joint limit m,n^ oo, -Ji ► p, we have the fol- 
lowing asymptotic expansion for the perturbed partition function Z m ^ n {ot, (3) . 
For all (a,/3) £l 2 , 



Z m , n (a,(3) = (-i)™» e -™/ 3 e- mnf 



/2^P( e -p7T)2 

x e ^ fe+ ^e-^ 2 /^ 2 )/ 2 (l + (l)), 

(k, 



where P(q) = YVk=ii^ ~ Q 2k ) an d f * s free energy per fundamental domain 
of equation (1). 

Remark 5. Recall that, by definition, the free energy per fundamental 
domain f is given by 

f = - lim — logZ mn (0,0). 

m,n— >oo fjin 

Hence, it is not surprising that f should govern the exponential growth rate 
mn of Z mn (a, /?). 
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Corollary 6. For all (a, 13) e M. 2 , 

lim F mn (a,f3) = F p (a,P). 

m,n— >oo r 
nl \/Zra— >p 

Proof. Combining Lemma 3 and Theorem 4, we have 

lim F mn (a,{3) = — V e <ok+M e -<# W 2 )/2 

m,n->oo ' Z p 

n/(y/3m)—*p ^ '' 



_L y e -K(ak+pe) e -n(k 2 / P + P e 2 )/2 

z 



E 

Jp (M)eZ 2 

= F p (a,P), 

where Z p = J2(k,e)ez 2 e _7r ( fc2//p+/rf2 )/ 2 . In the second equality, we have used 
the fact that the sum is symmetric in (k,£). □ 

Remark 7. Recall that F mn (a, (3) is the moment generating function of 
the random variable wind min , which is computed using the reference dimer 
configuration Mq. A natural question which arises is what happens when the 
reference dimer configuration is changed. Let F^f^^a, j3) be the moment gen- 
erating function of the winding number, computed using a generic reference 
matching M\. Then, looking at the proof of Lemma 3, we obtain 



F^ h n (a,(3) = 4^W)/(2m) B (-iVin(Mi)+JV II (M 1 ))/(2n) 



Z m , n (0,0) 

Using Theorem 4, we deduce that 

F^\(a,/3) = ANi(Mi)/(2m)-n/3 B (-N lll (M 1 )+N 11 (M 1 ))/(2n) 

x— V e 7r ( afc +^) e - 7r ( fc2 //' + ^ 2 )/ 2 (l + o(l)). 
Zp (fc/)ez 2 

Hence, Corollary 6 holds provided that the reference dimer configuration 
satisfies JVi(Mi) = ^qp and iV m (Mi) = N n (Mi). Observing that N h N n , 
Nui are always constrained to satisfy N\ + Nn + Nn = 2mn, this implies 
that JVi(Mi) = JVn(Mi) = JV ra (Mi) = 2ffi2. 

3. Proof of Theorem 4. The main ingredient in the proof of Theorem 
4 is an explicit expression for the perturbed partition function Z mn (a,/3) 
as a combination of Jacobi theta functions. This is given in Proposition 8, 
Section 3.2 below. The proof of Proposition 8 is postponed until Section 4. In 
Section 3.1, we recall the definition of Jacobi theta functions. Proposition 9 
of Section 3.3 gives a concise formula for the combination of theta functions 
involved in the expression of the perturbed partition function Z mn (a,f3). 
Section 3.4 consists of the proof of Theorem 4, using all of the above. 
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3.1. Jacobi theta functions. Recall the definition of Jacobi's four elliptic 
theta functions "&i(C, 9), i = 1, ... ,4, and their expressions in terms of infinite 
products: 



MC,q)= E (-l) fc -V2 e (2*+iKg(*+i/2)- 

k=— 00 

00 

= 2g 1 / 4 sin(C)P((/) 11(1 " 2<Aos(2C) + 
1=1 

00 

& 2 ((,q)= ]T e (2fc+i)* g (*+i/2)' 
fc=— 00 

00 

= 2gV4 cos (C)P(g) J] (1 + V cos(2C) + 



00 



*s(C,?)= E e 2fc V 2 =^)Il( 1 + 2 ^ +lcos ( 2 0+^ +2 ), 

fc=-oo 1=0 

00 00 

*4(C,?)= E (-l) fc e 2fc V 2 = J P(^)n( 1 -V m cos(2C)+^ +2 ), 

fc=-oo £=0 

where -P(<?) = n^Li(l ~~ Q 2k )- It is sometimes convenient to use the notation 
^(C|r) for ^(C,e^ r ). 

The proofs of these formulas, as well as other properties of Jacobi theta 
functions, can be found in [2]. 

3.2. Perturbed uniform partition function. Let us recall the definition of 
the perturbed partition function Z mi „(a,/3). It is the partition function of 
the dimer model on the graph H m>n , where edges are assigned perturbed 
uniform weights: edges of type I (resp., II, III) have weights a = e ~ an /( 2m ) 
(resp., ft- 1 = e-^/O), b). Recall, also, the notation A = e a7T , B = e^ and 
p = lim m)n _ >00 ^t|^- An explicit formula for the partition function Z mjn (a, ft) 
is given in [8]: 

(3) 

+ Z^(a,f3) + Z^(a,(3)), 



where 



z& (<*,&= n n 

z ™=(-iy w n ={-i) r i 
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(4) 



P(z, w) = det 



1 b 

I + - 

w 

az 



b + 



w 



1 b\f, w 
b + w)[ b+ T 



a 2 z ■ 



w b 2 
tr w 



a 2 z. 



Each of the four terms Zm$ (oc, f3) is a determinant of a Kasteleyn matrix, a 
relative of the adjacency matrix of the graph H mn , computed with discrete 
Fourier transforms, using the invariance of the graph under translations by 
the vectors x and y of Figure 1. Using Riemann sums and analyzing the 
behavior of P(z,w) in the neighborhood of its zeros, one deduces (see, e.g., 
Kasteleyn [8]) that the free energy per fundamental domain with weights 1 
on the edges is 



f : 



1 If 

lim logZ mjn (0,0) = — -r \ log 

i,n^co ran iir z Jj 



1 

w+- + 2 

w 



\ dz dw 
Z ~ 



z w 



I pz7t r" log(2(cos ^ +1) _ e ;^ d( ^ 

47T Z Jo Jo 



r 2ir r27r 



which is precisely equation (1). 

The following proposition gives the asymptotic expansion of the four 
terms Zm$(a,f3) involved in the explicit expression (3) of Z m ^ n (pt, j3). The 
proof is postponed until Section 4. 



Proposition 8. In the joint limit m,n — ► oo, 



p, we have the 



following asymptotic expansion for the four terms involved in the partition 
function Z m>n (a, 0) : 



-i\n 7 (00) 
" J. ) ^m,n I 



("I) 



-iyi 7 (01) 
7 (H) 



z. 



(10) 



a, f3) 
a, /3) 
a, (3) 
a, (3) , 



> = (-l) mn A- n / 3 e 



7T0 2 p/2 p Trp/6 

-mnf ° ° 



P{qf 



(-M(,q)MCq)) + o(i), 

MC,q)MC,q) + o(l), 
where £ — § (pa + i(3), q = e~ pn and f is given in equation (1). 



3.3. Recombining Jacobi theta functions. The following proposition gives 
an explicit expression for the combination of Jacobi theta functions involved 
in the perturbed partition function. 
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Proposition 9. 

VC = x + iy G C, Vr G M = {z G C; Imz > 0}, 



E**(Cb-)ft(cir) 



i=l 



/2 * e -i2x2/(^r)^ 

r 



1 

"27 



i?3 [w 



Proof. Let q = e l7TT . By definition of Jacobi theta functions, we have 

MC\t)M(\t) + M(\t)M(\t) 



k+i\A(2k+l)C p i(2€+l)C „(fc+l/2) 2 + (m/2) 2 



fc+A_j2(fc+m)z J2(k-l)iy (fc+l/2) 2 +(£+l/2) 2 



(fc. 

The general term of this sum is nonzero if and only if k ^ £mod2. In this 
case, we introduce the two integers u and v such that k + £ + 1 = 2u and 
k — £ = 2v + l, and rewrite #i(C|t)i?i(C|t) + i? 2 (C|t)^2(C|t) as 



(5) 



M(\r)MC\r) + MC\r)MC\r) 

_ 2 e i2u ( 2x ) e i2(2v+l)iy ^2« 2 +(2d+1) 2 /2 

(it.-y)eZ 2 

Similarly, we write 

^3(C|r)^3(C|r)+^4(C|T)^4(C|T) 

= J2 (l + (-l) k+e )e l2k <e m <q k2+i2 

(fc/)£Z 2 

_ (1 + (-_ 1 )fc+^ e i2(fc+Qx e i2(fc-f)iy g fc 2 +£ 2 _ 

(M)GZ 2 

In this case, the general term is nonzero when k = £mod2. Setting u = 

k—f 

and v = —y- , we get 



(6) 



M(\t)M(\t)+M(\t)M(\t) 

_2 e i2u(2x) e i2(2i))i?/^2« 2 +(2i)) 2 /2 

(u,w)ez 2 

Summing (5) and (6), and recombining the sum over v, we obtain 



E^(C|r)^(C|r) = 2 E e^^e^^'^ 2 

(u,v) 



t=l 
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2 g e 2iu(2x)^\l g u 2iv(iy)( q l/2y* 



2i9 3 (2x|2r)i? 3 Uv 



"27 )<M *l/ 



The last equality results from the Jacobi identity which describes the trans- 
formation of the function #3 under the modular group [2] : 



Vu £ C, Vex £ HI t? 3 ( u |ff) 
with u = 2x and a = It. □ 



Corollary 10. In the case where £ = ^ (pa + i(5) and r = ip (i.e., 
q = e-*P), 

4 it: 



i=l 



P 



lira 
'~2~ 



1 
2^ 



iir/3 



-na 2 p/2 



J2 A k B 



k jyi n (k 2 /p+p£ 2 )/2 



where A = e™ and B = e n/3 . 



3.4. Proof of Theorem 4- Using expression (3) for the perturbed parti- 
tion function Z m ^ n (a,(5) and Proposition 8, we know that in the joint limit 



m,n — > 00, 



\f?>m 



Pt 



Z m!n (a,(3) = (-l) mn A~ n ^e- mnf 



e ira 2 p/2 e iTp/6 



(7) 



2P(q) 2 



u=l 



Replacing the combination of Jacobi theta functions of (7) by the expression 
given in Corollary 10 yields 

Z mtn (a,P) = (-ir n A- n / 3 e- mnf - 



(8) 



!TpP( q y 

x Yl A k B"e^ k2 IP + P^I 2 {l + o{\)) 
(fc,f)ez 2 

which is precisely Theorem 4. 
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4. Asymptotic expansion of the perturbed partition function: Proof of 
Proposition 8. Let us quickly recall the content of Proposition 8. It gives 

explicit expressions for the four terms Zm°n{a, /3), Zm,ri(a,/3), Zm]n (ct, /?) , 
Zm°n (a, (3) involved in the expression (3) of the perturbed partition function 
Z mtTl (a,j3), as functions of Jacobi's four theta functions. Sections 4.1-4.4 

consist of preliminary computations for each of the four terms Zm$ (a, (3). 
Proposition 8 is then proved in Section 4.5 using the aforementioned compu- 
tations. Let us recall the notation used: a = e -W( 2? ™) ; b = e ^/3/( 2n ) ; A = e na , 
B = e *f , ,p = lim™,™ C = |(pa + i0) and q = e~^. 

4. 1 . Computations for Z^jn i a iP)- The beginning of this computation is 
inspired by [14]. Using the explicit expression (3) for z£$(a,P) yields 

z<$(a,0)= n n (S+( 2 -° 2 *)+£ 

2 ™=-i«i»=-i 



-1 w n =-l 



b 2 



w 



n n ^fS-^OfS-^ 



1 w 

-1 w n =— 1 



II II ^-^-b 2 n)(w-b 2 r 2 ), 

™=-l w n =-l 



where 




a 2 z 



ri ,r 2 = -l + — ±\ 1-— -l. 



Since n and r 2 depend on a and z through the product a 2 z, we define <p G 
as a 2 z = e*^ and 




1_ 

Recall that a = e _7rQ /( 2m ) anc j that z is an mth root of —1. Hence, for large 
m and n, we are concerned with values of 4> whose imaginary part is close 
to and whose real part is in [— tt, it] . 

Choice of r\ , r 2 and their domain of definition. In the formula for n ((j>) , r 2 (0) , 
there is an ambiguity in the determination of the square root, which we 
now clarify. The functions n((/>), r 2 (<^>) are analytic on a domain D, pro- 
vided the square root does not vanish in D. The square root vanishes when 
(Re(0),Im(</>)) = (0[2vr],log(l/4)), so let us choose 

D = {(f> e C|-7r < Re(0) < vr, -c < Im(</>) < c}, 
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Fig. 7. Plot of the two functions ri(<j>) (plain line) and r2(</>) (dashed line) for 
4> £ [— 7r,7r], and their position with respect to the unit circle. Note that when a is not 
and m is large enough, ri(i—) and T2{i — ) are still on the unit circle (see Lemma 11 
below). 

where c is some positive constant smaller than log(4). Now, observe that 
r i(<l ) ) r 2(4 > ) = 1 so that one of r\, ri has nonnegative imaginary part. More- 
over, n satisfies the polynomial equation 

(10) n + — = e** - 2. 
Taking imaginary parts on both sides yields 

(11) MriW))(l-u^)=e- Ini Wsin(Re(0)). 

Hence, Im(n (</>)) ^ whenever Re(0) ^ {— n, 0, ir}. When Re(^) = 0, we look 
at equation (9) and deduce that the term in the square root is positive when 
| lm(0)| < c so that Im(n(^>)) 7^ in D. Using the analyticity of r\ and 
in D, we conclude that one of r±, ri has imaginary part which is positive 
in D. By convention, we choose ri((f>) such that Im(ri ((/>)) > in D. The 
values for n, ri on the boundary of the domain D are chosen by continuity. 
Figure 7 shows a plot of r2(4>) when 4> is real. 

Let us first compute the product over w in (8). Since w is an nth root of 
— 1, we have 

VAgC Yl (A - w) = X n + 1. 

In particular, 

IJ ^=(-l) n , II (#r j (<l>)-w) = l + Br?(<l>), j = l,2, 
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so that 

(12) Z^(a,P) = (-l) mn B- m J] (1 + U7?(0)(l + Utf(fl). 

j.m = _l 

Ideas governing the next steps of the computation. There is a change of 
behavior in the product (12) when |rj(</>)| is smaller, greater or close to 1. We 

start by factoring out rj(<j)) when > 1. This part of Zm,n(ct,/3) grows 

exponentially, so we expect it to involve the free energy per fundamental 
domain f. We postpone its study until Section 4.5. Then, in Lemma 12, 
we deal with the part of Zm}n(ct,(3) where |r.,(</>)| < 1. This part makes no 
contribution to Zm,n(a,/?)- Finally, in Lemma 13, we deal with the part of 
Zm^n(oi, 0) where |rj(0)| is close to 1. This part involves the Jacobi theta 
functions giving the discrete Gaussian term of Zm,n (a, /?). In Sections 4.2, 
4.3 and 4.4, we study, in a concise way, the computations for Zm°n(a, (3), 
^m,n(a,/3) and Z^n°n{a, (3). 

In order to proceed with these steps, we need: 

1. a characterization of > 1, |r,-(0)| < 1 and = 1; 

2. an expansion of rj{4>) when cj> is close to (for Lemmas 12, 13); 

3. to know that is increasing as a function of Re(^>) when — ir < 
Re(^>) < and Im(<^>) is small (for Lemma 12). 

All of this is provided by the following lemma. 

Lemma 11. On the domain D, the functions n andr 2 have the following 
properties: 

1. - // < Re(0) < tt, then |n(0)| > 1, \r 2 {4>)\ < 1. 

- // -7T < Re(0) < 0, then |n(0)| < I, \r 2 (i)\ > 1. 

- J/Re(0) = O, then |n(^)| = \r 2 (<f>)\ = 1. 

2. Their expansions in the vicinity of(p = are given by 

(13) n (0) = e i27r / 3 (l + A + 0(02)^ = e an/3 e */V*H>tf) > 

(14) r 2 (<^) = e" 12 ^ 3 (l " ^ + 0(<A 2 )) = e - i2 -/3 e -0/V3+O(^)_ 

3. When — ir < Re(</>) < 0, |ri(<^)| is increasing as a function o/Re(</>). 

Proof. 1. Since r2(0)?r(^>) = 1, it suffices to prove the properties for 
n (</>)■ Recall that, by convention, lm(n(0)) > in D. Looking at equation 
(11), we deduce that when sin(Re(0)) > or, equivalently, Re(0) G (0,7r), 
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we have |n(<^)| > 1. We further deduce that when sin(Re(0)) < or, equiva- 
lently, Re(0) G (— ir, 0), we have < 1- The explicit evaluation of ri(±7r) 

shows that the inequalities are still strict when Re(^>) = ±7r. 

We thank the referee for the idea behind the above argument. 

Let us now consider the case where Re(^) = 0. Since <j) G D, we have 
lm(0) >log(l/4) and 



with no ambiguity in the square root. An explicit computation yields |ri($>)| = 
1. 

2. The second part of the lemma is obtained via an explicit Taylor ex- 
pansion of ri((p),r2(4>) in a neighborhood 0, which can be performed since 
these functions are analytic. 

3. We need to prove that ^r^I) * s n o nn egative whenever — tt < Re(0) < 
0. The derivative of |ri(</>)| 2 with respect to Re(</>) is 



From this, we deduce, by differentiating with respect to Re(</>), that 




so that, in this case, r\{(j)) is given by 






dn(4>) 
<9Re(0) 

9|n(0)| 2 



ie^ri(0) 




2ri(0) + 2- e^' 
/ ie^\n(0)\ 



2 



Hence 




= sign(lm(e-^(ri(0) + l))). 
Rewriting the quadratic equation satisfied by r\{<f>) as 

(ri(0) + l) 2 = e ^ri(0), 
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we see that e l( ^(ri(<j)) + 1) equals pj^r^ , whose imaginary part has the 
same sign as Im(ri (</>)), which is nonnegative by definition. □ 

Since z is an mth root of —1, we can write z = e *>r(?3+ 1 )/™ for j G {— [^J > • • • > 
L^J} so that = ^(*« + 2 i + !)• Using Lemma 11, we know that when 
j > 0, |?r(0)| > 1, and when j < —1, | T2 (0) | > 1- Factoring those terms out in 
(12) and using an algebraic manipulation to produce a product over positive 
j's only, we obtain the following expression for Zm}n(ct, (3): 



(-1)">"B-'" L '"jJ ( 1 + Br n (*(><* + 2 J + !' 



i=-Lm/2j 



??? 



x(l + Sr^ 7r(?a + 2j + 1 ^ 



|(m-l)/2| , ,. 

/ 7r(za + 2j + i; 



(-l) mn II r i 
V i=o 



■/?? 



x 



t^- 1 (2j + l)) 



n ^ 

\ j=0 
L(m-l)/2j 



vr(ia + 2j + 1) 



x 

3=0 



x ^1 + Br£ 

|m/2j-l 



■/?? 



X 

3=0 



l + B ^ n (^a-{2 3 + l)) 



m 



Lemma 12. In the joint limit m,n^ oo, -J== > p, we have 



j=[V /4 J 



m 
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L(m-l)/2j 



j=LmV4j 



n u+ s 



ir(ia + 2j + 1) 



■/?? 



|m/2j-l 

n u+s 

j=LmV4J 



1 -n/^ 7r ( i «-( 2 i + 1 )) 
7*< 



■/?? 



L TH<«^)H + » (1 ,, 

j=L m 7 J 



l + o(l). 



Proof. We give the proof for the first product on the second line. The 
proof for the other terms is obtained using the same arguments. 
Taking the logarithm of the left-hand side, we get 



[m/2j-l , 
£ log! 1 + Brl 

j=Lm 1/4 J 

[m/2\-l 



n fir(ia-(2j + l)) 



in 



< E 

[m/2\-l 

< E cB 



log 1 + Br? 



B /7r(*a-(2j + l)) 



??? 



vr(ia - (2j + 1)) 



for some constant c> 



< mcfi max 

\m 1 / A \<i<\m/2\-l 



n f7r(ia-(2j + l)) 



m 



< mcB 

< mcBe 
= o(l). 



ix(ia - 2m 1 / 4 ) 
m 

-27rpm 1 / 4 (l+o(l)) 



1 



(Lemma 11, part 3) 
(Lemma 11, part 2) 



□ 



Lemma 13. In the joint limit m,n^> oo, — ► p, we have 

'7r(*a-(2j' + l)) 



Lm n J "Vi +5 -vr 



3=0 



m 



1 + Br^ 



m 



+ o(l), 



m / / V V m 
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n 

J'=0 

^3(C,g) , m 

= ^T + o(1) ' 

w/iere C = § (pa + «/3), q = e~ p7T and P(q) = n£?(l ~ ? 2fc )- 

Proof. The arguments used are similar to those of Lemma 12. First, 
note that 

oo 

MC, q) = P(q) 11(1 + e 2 V i+1 )(l + e"^ 1 ) 
j=0 

= lim P(g) TT (l + e 2 V J ' +1 )(l + e _2 V J ' +1 )- 
j=0 

Let us prove the following (the other three cases are handled similarly): 



n (i+s-v 

(15) 

\m L / A \-l 

= II (l + e 2 V +1 )+o(l). 

j=0 

A generic term of the right-hand side of (15) can be written as 

(1 + e 2 V i+1 ) = (1 + e~^e i7Tap e- np{2j+1) ). 
Note that such a term can be equal to 0, whenever a = 2 +1 for some integer 

k and f3 = —p(2j + 1) for some j G {0, . . . , [wi 1 / 4 ] — 1}. So, the product of 
the right-hand side of (15) might vanish. Observe that when this is the case, 
there is only one term in the product which is actually 0; let us denote by j 
the corresponding index. 

We first consider the part of the product corresponding to nonzero terms 
on the right-hand side of (15). That is, depending on the values of a and /3, 
we either consider the whole product or the product over all indices except 
j. In order to keep the notation simple, it is implicit that j is omitted 
whenever it occurs. Dividing the left-hand side by the right-hand side in 
(15) and taking logarithms yields 



n(ia - (2j + 1)) 



E Jos 

3=0 



1 + fl-Vg - (2j + l))/m) 

1 + B- 1 e i7Ta Pq 2 ^ +1 
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|_ m V4j_l 



< 



E 

3=0 



log 



1 + B~ l r^ n (-K{ia - (2j + l))/m) 



l + S- 1 e i7ra ' , g 2 -5 , ' +1 



Observe that the two quantities B^e^q 2 ^ 1 and fl-i r -"( 7r ( iQ -(y+ 1 )) ) 
are bounded, and bounded away from —1, uniformly in j = 0, . . . , |_n2 1//4 J , 
j ^ j, for m and n large enough. We therefore get the following bounds, 
where c\,C2 are positive constants: 



LmVj-i 

E ^ 

j=0 



1 + B- l r2 n {n{ia - (2j + l))/m) 



[mV^J-l 

< E ci-b- 1 

i=o 



1 + B-l e «™*Pg2j+l 
n fia-(2j + l) 



■m 



g inap 2j + l 



L m V4j_l 

< ^ Cl . B -l| e 7r(ia-(2j+l))n(l+o(l))/(v / 3m) 
i=0 



_ e 7r(ia-(2j+l))pi 



Im 1 / 4 ] 

<c 2 ^ |ia-(2j + l)|e"^ +1 )^ + °( 1 » 
i=o 



by (13) 
n(l + o(l)) 



because the series with general term |ia — (2j + l)|e' 



vom 



-7r(2j+l)p 



0(1) 



is convergent 



and thus the upper bound goes to zero when goes to p. 

Let us now consider the case when a = 2*±1 f or SO me integer A; and /3 = 

—p(2j + 1) for some j £ {0, . . . , [m 1 / 4 ] — 1}. For all terms involving indices 
j j, the above estimate holds. When j = j, the corresponding term on the 
right-hand side of (15) is 0; let us show that the corresponding term on the 
left-hand side is o(l): 

" V m 

= |1 + 5-i e ^(ia-(2j+l))«(l+o(l))/(v / 3m)| by ^ 

= 11 + 5- 1 e 7r ( iQ -(2i+i))p[ e 'r(^-(2i+i))(n(i+o(i))/(v'3m)-p) 

ra(l + o(l)) 



l+7r(io- (2j + l)) 



■m 



P)+o(l) 



ir(ia - (2j + 1)) 



n(l + o(l)) 



y/3- 



m 



P)+o(l) 



o(l). 



LOOP STATISTICS IN THE TOROIDAL HONEYCOMB DIMER MODEL 23 

In the two last lines, we have used the assumption 1 -\- B~ 1 e n ^ ia ~^ +1 ^ p = 0. 
Combining this with the estimate when j ^ j yields the lemma. □ 

Corollary 14. In the joint limit m, n — > oo, -7=^ — > p, we have 



ALJa)ALn(a) P{q) 2 



where 



A k»=[ n r \- — - — ■ 



in 



4.2. Computations for Zm°n{ot,(3). The computations for Zm°n{a,(3) go 
through in a similar way. Let us just stress where the differences occur. Using 
the explicit expression of (3) for zgJl(a,/3) yields 

zi$w)= n n ^-(^-^i)(^-fe 2 r 2 ). 

This time, w is an rath root of 1 (instead of —1), so 

VAGC Yl (X — w) = \ n - 1. 

10" = 1 

In particular, n«j n =i w = ( — l) n_1 and, when performing the product over 
w, we obtain 

x l( "rf J f 1 - Br?( M"»+2 J + im f 1 _ Br? M» + + 1: 



j=-[m/2j x 

/L(m-1)/2J 
(_!)mn JJ ri 

V j=0 
/Lm/2J-1 

X [] r 2 

\ j=0 

L(m-l)/2j 

X 

3=0 



m / J \ V 771 



ir(ia + 2j + 1] 



in 



Tr(ia - (2j + 1)) 



■in 



ii-B~\r nf<ia+23+l] 



in 
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'n(ia + 2j + 1) 



x I 1 - Br r i 

|m/2j-l 



/7? 



X 

3=0 



n H<^^)) 

vr(ia - (2j + 1)) 



x ( 1 - B- X r 2 n 



m 



The rest of the computation goes through in the same way, except that the 
+'s of Zm : n are replaced by the — 's of Zrn^n- As a consequence, the analog 
of Lemma 13 involves the fourth Jacobi theta function ^(Cw) instead of 
the third one. We summarize the expression for Zm^n in the following result. 

Corollary 15. In the joint limit m,n^ oo, -7=^ — ► p, 

(-ir n z^(a,(3) _ M(,q)MC,g) o(1) 

A^,Ja)A^ n (a) P(q) 2 

4.3. Computations for Zm)n (a, (3) . The computations for Zm]n (a, (3) are 

slightly different. When performing the product over w for Zm)n(ct, j3), we 
get 

L(m-1)/2J , r , „, v 

Z^(a,P) = (-l)™B-™ J] ( 1 + gr ?( m ), 



j=-[m/2\ 



m 

In order to obtain an expression similar to what we had above for Zm}n and 

Zm°n, we isolate the term j = and factor r^(^-(ia + 2j)) when its modulus 
is greater than 1. This yields 

(16) 



lair \ 
m J 



'L(m-1)/2J ,. \n/|m/2j ,. , 



L(m-l)/2j 



x n ( 1 + 5-Vr" 



7r(«a + 2j) 
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' i\(ia + 2j) ' 



x ^1 + 

\m/2\ 



m 



(17) x J] [ 1 + Br-J <** ~ 2j) )) ( 1 + B- V^< fa " 2j) 



m \ \ m 



Let us rewrite the first two terms in brackets, using the fact that r±(^-) 
„ 1 „ / iair \ \ ( „ „ / ia7r 



(18) 



S -l/2 rr /^^ +jB l/2 r n/2/^ry 



: 4 cosh (log (l? 1 ^ (^) ) ) cosh (log (* V2 r »/2 (^ 



and note that 



hn V 4 cosh (log (^ 1/2 ^ /2 (^))) cosh (log (bV\ 2 ^(!^ 



n,m— >oo 
n/(V3m)->p 



.^-z^Wvr^ + z/?) 



2 

(2cos(C))(2cos(C)) 



For the other terms of (16), we use the same arguments as for Zm}n(a,f3) 
and obtain: 

Corollary 16. In the joint limit m,n^ oo, -7=^ — > p, 

+ o(l), 



r m,n(a) r m,n(") P (<?) 

/l(m-l)/2l , ,. 

r m,n(a)= 11 n( 

V j=l V 



???. 



V j=i 



LW2J / /• -\ 
/ 7r(iQ! — 2j) 
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4.4. Computations for Zm°n (a, (3) . Computations are similar to those for 
Zm}n(oi, (3). The analog of (18) is, in this case, 



1 ~n(^\\( 



B - r ^)){ 1 - Br ^-J 



= 4 sinh (log (j3V2 r »/2 (^J) si nh (log (b 1 ' 

which, in the scaling limit, converges to 

4s . n jMia^\fn(-ia P + (3) 



\ m J 



2 J V 2 
= lA ^-^±) S m(^^±M>) = (2 8 in(C))(2 S in( C )). 

As in the case of Zm,i! , the product over the roots of unity introduces a factor 
(— l) m ( n_1 ). An extra factor (— l) m_1 appears since we factor —rk{ ' K<y%a ^ 2 ^ ) 
when its modulus is greater than 1, something which occurs exactly L^^^^J + 
[?rj = m — 1 times. Thus, we have the following result: 



Corollary 17. In the joint limit m,n^ oo, — > p 



+ o(l). 



4.5. Investigation of K l m n (a) and T % mn (a), and proof of Proposition 8. 
As mentioned in "Ideas governing the next steps of the computations" in 
Section 4.1, we deal here with that part of each of Zmjn(ot,0), zi^a, (3), 

Zm}n(a,f3) and Zm°n.(a,{3) which grows exponentially in ran. Recall that 
we expect this growth rate to be driven by the free energy per fundamental 
domain f of equation (1). Proposition 8 below gives precise statements about 
this by exhibiting the relation between A^ n (a)A^ n (a), n (a)r^ n {a) 
and A^ n (0)A^ n (0), and giving the asymptotic behavior of A^ n (0)A^ n (0). 

Note that the proof of Proposition 8 is completed by combining Corollaries 
14, 15, 16, 17 and Proposition 18. 

Proposition 18. In the joint limit m,n — > oo, — >■ p, 

1 An/3 ^rn,n( a )^-m,n( a ) ira 2 p/2 1 1 _i_ 1 \\ , 

l - ■ A^JOJA^JO) ~ e -U + oUJb 

9 An/3 (- 1 )" r m,n( a ) r m.n( a ) 1/2 c wa 2 p/2 n-L^n^- 

2 ' A ' Ai,,„(0)A;L n (0) ~ 9 ' 6 + 

3. A^ n (0)A^ n (0) = e^ ■ e~ mnf • (1 + o(l)). 
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Proof. Recall the definition of the functions r±((f)) and T2((j>) given in 
(9). Recall, also, that there is an ambiguity in the choice of the square root, 
which we had clarified on the domain D = {<p G C|— ir < Re((/>) < it, —c < 
lm(0) < c}, where c is some positive constant smaller than log(4). 

For the purpose of proving Proposition 18, we need to extend this clar- 
ification to a domain D' containing — ir < He(4>) < 3tt and small values of 
Im((j)). On this new domain, we also want to give a meaning to log(ri (</>)). 
Hence, we need to show that we can associate with r\((f)) a well-defined 
argument in (0,27r). 

Recall that ri((f>) is analytic provided that the term in the square root of 
(9) does not vanish. This only happens when (Re(4>), Im(^)) = (0[27r], — log(4)), 
so the function r\ is analytic on 

D' = {(j) e C|— 7r < Re(» < 3vr, -d < lm(0) < d}, 

where d is some positive constant smaller than log(4). 

Equating the real parts of the two sides of equation (10) yields 

Re(n(0)) (l + r^yp) = e- ImW cos(Re(<A)) - 2 

so that Re(n(^>)) < on D'. As a consequence, one can associate, without 
ambiguity, an argument in (0, 2ir) [even in (§,%-)] with r±((j)). We deduce 
that the function log(ri) is analytic on D' . 

We need one more property of r\ : looking at the term in the square root 
of (9), we know that r\{(f) + 2ir) = r2((f)) = r /^\ ■ Therefore, 

(19) arg(n(0 + 2vr)) = 2ir - arg(n(0)). 

Proof of 1. Let us first consider the log of the ratio ^T'" T' n [^} ■ Recall 
that r2((p) = r\((j) + 2tt) so that one can write 

/|(m-l)/2| , „ . 

A' (a)Al n (a) =- 1 




It follows that 



, KnMAlM xrVi /7r(za + 2j + l)\ /vr(2j + l)^ 

A m,n(°) A m,n(0) V m 
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m— 1 

3=0 

where /i(0) = logn(0), ^ = and e Q = EES. We know that /i 

analytic in Z)', so we can perform a Taylor expansion. This yields 



, A m,n( Q! ) A m,n( Q! ) 

m-1 / 2 

For the first term, we have 

j=0 J0 3=0 J <t>j-n/™ 

m ~ 1 r<t>j+n/m 



/ 1 (2tt)-/ 1 (0)- ^ / " /f(^.)(0-^)# + O — 

„_n J<t>i—Tr/m \7Jl 



/i(27r)-/i(0) + of 



3=0 -"Pi-^h 

m 2 



For the second term, 

2tt ^ / >27r / 1 \ / 1 

X "* nil / , \ I nil / ,1 , . I L \ nl ,„ \ nl I „\ ^ / X 



3=0 



/-27T / 1 \ 

E/iU)= / fittW + O (-J =/l(27T)-/((0) + O 
;_n ■'O \Ul J 



Recombining the different terms gives 

-^ (/!(2 .,- /!( o,) + o(^), 

where / 1 (2tt) = ^f, A(0) = *f , /((2tt) = /((0) = ^. 

Hence, 

log tt — / x ■ o /x = log A H h G> — ^ 

which proves part 1. 
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Proof of 2. To prove part 2, it suffices to show that lim^m—^ (— l) n A^ n (a) 

X ^m,n( a )/ T kn( a ) T2 m,n( a ) = = ^ ■ We Write T m,n{a)T 2 m n {a) as a 

unique product: 

4(m-l)/2l , Im/2J ... „ _ . N 

/ Ti{ia + 2j)\ L J 4 J /7r(?a + 2m-2j) 



n r, n 



3=1 3=1 



''1 



/ 7r(ia + 2j)\\ /"j-r 1 / fit{ia. + 2j) 



(20) =n»FP = n 



n 



,j=l / \j=l 

'7r(ia + 2)\ {Tr(ia + 2(m - 1)) 



m J \ m 



*nf(«M«^" 

The first two terms of the last line consist of the first and last terms of the 
product (20), while the product over j contains the terms of (20) grouped 
in pairs, starting from the second one. 

From equation (19), we know that ri(<j>)ri((j) + 2ir) = 1 and arg(ri(</>)ri(0 + 
2-71")) = 2-7T. In particular, we deduce that 



m J \ m 



Therefore, 



Let us consider the log of the ratio tztw^^II^^^TTa • K can be rewritten 
as 

log 



(-i)™ri lin ( a )r^ n (a) 

_ n .^i/ /i /^^x_ 2/iM ^ + 2 , + i; 



2\ \ m 1 \ m 

3=0 



ir(ia + 2j + 2) 
m 
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Since the function fx is analytic in D' , a Taylor expansion yields 



log 



Am \m J 



3=0 



4m Jo 



2tt 



n 



m- 



n 



Hence, 



log 



^(/((2.)-/((0)) + o(^). 



(- 1 ) nr m,n(«) r m,n(a) 4 ™ 

n7r 



2^: 



m 



n 



m- 



+ ^ ='- + ^ , 



n 



m- 



which proves part 2. 

Proof of 3. Recalling the definition of the free energy per fundamental 
domain f given in equation (1), let us integrate it explicitly over ip: 



f : 



I r2n r2n 
'4^2 



log(2(cosV + l)-e^)#(^ 

— r / / log \2(cos ip + 1) — e l ^\d(p dip (by symmetry) 
Jo Jo 



J 

2-7T t"K 



1 

"2^2 

VT Jo 
1 



7T r27Y 



.70 



7T r27Y 



log\e^ -rx{cf))\ dip d<p+ / / log - r 2 (^)| ## 



jo 



log|ri(0)| d<p 



2?r 



log rx (<f>)d<j> + / logr 2 ( 



2- 



27T Jo 

where, in the fourth line, we have used the identity 
1 r 27T 



2vr 



log if \t\ > \s\ 

log |s|, if |s| > |t| 



log|t + se^|dV : 
Let us now consider the logarithm of A^ n n (0)A 2 re n (0): 
logA^ B (0)A^ n (0) 



m— 1 
J=0 
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mn 
mn 



2k 



m ~ 1 r<f>j+w/m 



j=0 J<t>j-Tt/m 



2tt 



-^(/i(2")-/;(o))+o(4 

12m \m J 
by the Euler-McLaurin formula. Hence, we deduce that 

TV i TL \ 

log ALn(0)A' (0) = -mn . f + _ p + O _ . 
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